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1. Introduction 

In recent years we have witnessed an intensive research activity aiming at improving our un- 
derstanding of the behaviour of Yang-Mills theories in the low energy (infrared) regime. This is 
a rather complicated issue, closely related to the confinement of quarks and gluons. We consider 
pure Euclidean SU (N) Yang-Mills theories with action 

Sym = ^ j d A x F^ V F^ V , (1.1) 
where A^, a = 1, ...,N 2 — 1 is the gauge boson field, with associated field strength 

F% = d,A a v - d v Al +gf bc A b ^y . (1.2) 
The theory is invariant with respect to the local gauge transformations 

8A a ll =Dfco h , (1.3) 

with 

Df = d^ ah -gf ahc A^, (1.4) 
denoting the adjoint covariant derivative. 

Thanks to the asymptotic freedom, the coupling constant of the theory turns out to be small at very 
high energies, where perturbation theory then becomes reliable [[j], |2p. However, in the infrared 
region, the coupling constant grows and nonperturbative effects have to be taken into account. 

The introduction of condensates, i.e. the (integrated) vacuum expectation value of certain oper- 
ators, allows one to parametrize certain nonperturbative effects arising from the infrared sector of 
e.g. the theory described by (|Q|). Via the Operator Product Expansion (OPE) (viz. short distance 
expansion), which is applicable to local operators, one can relate these condensates to power cor- 
rections which give nonperturbative information in addition to the perturbatively calculable results. 

Possible causes of nonperturbative effects are given by Gribov ambiguities which affect the Faddeev- 
Popov quantization procedure and hence the propagators pp, the existence of (topologically) non- 
trivial field configurations like instantons H], etc. These effects are not necessarily unrelated, as 
e.g. instantons can contribute to the condensate (a 4 .F^ v ) [Q], the Faddeev-Popov operator has zero 
modes in an instanton background etc. 

Although propagators are not gauge invariant quantities, they are, to some extent, the simplest 
Green functions which can be studied from the analytic point of view. As far as four-dimensional 
Yang-Mills theory is considered, this task appears to be very difficult for more complicated Green 
functions. 

During the last decade, there has been an intensive activity from the lattice community in the study 
of the gluon and ghost propagators. We can thus compare, at least qualitatively, our theoretical 
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predictions with the available lattice data. 

So far, a certain number of gauges have been considered extensively from the theoretical as well 
as from the lattice point of view. This is the case for the Landau, Coulomb and maximal Abelian 
gauges. Throughout this paper, we shall mainly be interested in the Landau gauge 

<9M*=0. (1.5) 

In particular, we shall focus on the effect of the dimension two condensate (A^ in ) on the infrared 
behaviour of the gluon propagator, in combination with the effects arising from a treatment of the 
Gribov problem. 

2. The dimension two condensate 

The possible existence and relevance of a gauge condensate of dimension two has been advo- 
cated a few years ago by [§,0]. The idea was put forward that the gauge invariant, dimension two 
operator A^ in , obtained by minimizing A 2 along its gauge orbit, 



— fd 4 x(A u u ) 2 
/T J v V 



UeSU(N) VT . 

A u = UA^ + Ud^ , (2.1) 



might condense, i.e. 



<AL>/0 (2.2) 
This nonvanishing condensate might be important for several reasons: 

• (A 2 ^ ) could serve as an order parameter for the condensation of monopoles, relevant for the 
dual Meissner effect in the dual superconductivity picture of color confinement [||, [7|] 

• it could account for certain power corrections in which have been reported in the study of 
two and three point correlation functions, see e.g. [§, 0]. 



it could give rise to a dynamical gluon mass m g via the OPE [10], which could be relevant for 



the dual Meissner effect [O] as well as for obtaining analytic estimates of glueball spectra 



[12]. 



3. A closer look at Ai 



mm 



The relevance of the operator A^ in for Yang-Mills gauge theories is known since many years. It 
plays a key role in the study of the Gribov copies and of the geometrical and topological properties 



of the space of gauge orbits [13, 14, |15j ]. It can be expressed as an infinite series of nonlocal terms, 
see [pJ|, namely 



^min I d X 



+ 0(A 4 ). (3.1) 
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From this expression, one can also check explicitly order by order in the coupling constant that 
A^ in is gauge invariant. However, it is immediately clear that A^ in is a highly nonlocal expression, 
a fact making it very difficult to handle. In general, it also falls beyond the applicability of the OPE, 
which refers to local operators. 

So far, the only possibility to study the operator A^ in by analytical tools relied on choosing the 
Landau gauge (|L~5|). Due to the transversality condition, all nonlocal terms drop out, so that A^ in 
reduces to the local operator A 2 , 

A min = A in the Landau gauge (3.2) 

3.1 Effective potential for the local composite operator 

The operator A 2 is an example of a local composite operator (LCO). In JT7|], a method was 
developed to construct a sensible effective potential for such operators. This so-called LCO method 
was later applied to the case of A 2 in the Landau gauge [|18|]. We couple the operator A 2 to the Yang- 
Mills action by means of a source J, 

Sj = S YM + J d 4 x Q/A£A£ - h^A • (3-3) 

The last term, quadratic in the source /, is necessary to kill the divergences in vacuum correlators 
like (A 2 (x)A (y)) for x — > y, or equivalently in the generating functional W(J), defined as 

e -W(J) = J [ fi ekis] e -^ 4j * . (3.4) 

The presence of the LCO parameter £ ensures a homogenous renormalization group equation for 
W(J). Its arbitrariness can be overcome by making it a function £(g 2 ) of the strong coupling con- 
stant g 2 , allowing one to fix £(g 2 ) order by order in perturbation theory in accordance with the 
renormalization group equation |l8| ]. 

In order to recover an energy interpretation, the term oc J 2 can be removed by employing a Hubbard- 
Stratonovich transformation 

1 = J oe~^^ + ? Al ~^ , (3.5) 

leading to the action 

S = Sym + So , 

s ° = S dAx {m + m goA2+ k {A2)2 )- (3 - 6) 

A key ingredient for the LCO method is the renormalizability of the operator A 2 . It was proven in 
Jl9| ] that A 2 is renormalizable to all orders of perturbation theory, making use of the Ward identities 
in the presence of the operator A 2 . In addition, an interesting identity was proven concerning the 
anomalous dimension y A i of the operator A 2 , first noticed in [|20|]. It can be shown that y A i can be 



4 



Infrared behavior of the gluon and ghost propagators in Yang-Mills theories 



Silvio Paolo Sorella 



expressed as a linear combination of the gauge beta function /3 and of the anomalous dimension Ya 



of the gauge field A, according to the relationship Q19[ ] 



y Al (fl) = -P + 7 A ( fl )j , <* = j^- (3.7) 

Starting from ( ^ ) it is possible to calculate the effective potential V(o). The correspondence 
(a) = — g(A 2 ) consequently provides evidence for a nonvanishing dimension two gluon conden- 



sate using an effective potential approach, if (a) ^ 0. It is clear from ( pffy that (a) ^ also induces 
an effective gluon mass. V(o) was calculated to two loop order in [18, 21], and a nonvanishing 
condensate is favoured as it lowers the vacuum energy. The ensuing effective gluon mass was found 
to be a few hundred MeV. 

3.2 Restriction to the Gribov horizon 

In the Landau gauge, it is necessary to restrict the domain of integration in the Feynman path 
integral at least to the so-called Gribov region £2, whose boundary dQ. is the first Gribov horizon, 
where the first vanishing eigenvalue of the Faddeev-Popov operator, 

Jt ab = -a M {d^8 ab + gf ac %) , (3-8) 
appears. This restriction is necessary due to the existence of the Gribov copies, which implies that 



the Landau condition (1.5) does not uniquely fix the gauge 



It has been discussed in [|22|, [23J how this restriction can be accomplished at the Lagrangian level. 



More precisely, the starting Yang-Mills measure in the Landau gauge is given by 

dfL y = DA8{d pL A a ll )fet{J?)e-( s ™ + ^ H ) , (3.9) 

where 

H = J d 4 xh{x) =g 2 J tfxf^A* (^- l ) ad f dec A e ^ , (3.10) 

is the so-called horizon function, which implements the restriction to the Gribov region Q.. Notice 
that H is nonlocal. The parameter 7, known as the Gribov parameter, has the dimension of a mass 
and is not free, being determined by the horizon condition 

(h(x))=4(N 2 -l) , (3.11) 

where the expectation value (h(x)) has to be evaluated with the measure d\iy. To the first order, the 



horizon condition (3.11) reads, in d dimensions, 



_N(d-l) 2 r d d k 1 

4 8 J (^MVgY ' 

This equation coincides with the original gap equation derived by Gribov for the parameter 7 [0]. 
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Albeit nonlocal, the horizon function H can be localized through the introduction of a suitable 
set of additional fields. The final action reads 

S = So-fg J d d x^ ahc A^ + f bc A a ^-d(N 2 -\)r 4 ) , 
So = S YM + J d d x (b a d^+^ (D^c) a ) 

+ J d d x(q>™d v (d v <p«; + gf abm A h v <p™) -T5 a ;d v (d v CD a ; + gf abm A b v co™) 

-g(d v W%)f bm (D v c) b <p™) . (3.13) 



The fields (lp a £ , (p'jfj are a pair of complex conjugate bosonic fields. Similarly, the fields (ooff, oo'jf) 
are anticommuting. The horizon condition is equivalent with the demand that the quantum effective 
action F obeys 

dT 

df = ° (3 - 14) 



As shown in p2 , 23[ [24 1, the resulting local action turns out to be renormalizable to all orders of 



perturbation theory. Remarkably, we have been able to prove that this feature is preserved when 



the local operator A 2 is coupled to the Zwanziger action [26]. This allows for a simultaneous study 



of the effects of the Gribov parameter and condensate (A 2 ). 

A main consequence of the restriction of the domain of integration to the Gribov region is the 
fact that the ghost propagator gets enhanced in the infrared region. Using the gap equation arising 
from the horizon condition, one finds that [|3[ ^2, 23|, 25, 26, ^7|] 



c a c h 



8 ab 

— for P 2 
P 







(3.15) 



This enhancement remains valid in the presence of (A 1 ) [25, Em. The infrared enhancement of the 



ghost propagator in the Landau gauge has also been observed from lattice simulations [29, 28] or 



solutions of the Schwinger-Dyson equations [30, 31, 32] 



The Gribov restriction and (A 2 ) also affect the gluon propagator in a nontrivial fashion, more 



precisely one finds [25, 26] 



A a ,A b ) 



tab 



PiiPv 



p 4 + m 2 p 2 + 2g 2 Nf- 



(3.16) 



It is worth noticing that a propagator like (3.16) is not new, as it has been considered already some 



time ago by Stingl in order to solve the Schwinger-Dyson equations [p3|. It induces an infrared 
suppressed gluon propagator, a fact in qualitative agreement with lattice [ p9| ] and Schwinger-Dyson 
results [Q |3~I, 32fl. Let us also mention that ( |3.16| ) violates spectral positivity, giving an indication 
that the gauge bosons are unphysical particles J26| ] 

4. Beyond the Landau gauge 

It is unclear what the role of A 2 - might be in other gauges. This is a very difficult question, 



without any answer at the moment. Due to the severe nonlocality in the expression (3.1), it seems 
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Gauge 


LOperator 


linear covariant 
Curci-Ferrari 
maximal Abelian 


1 A a A a 

2 H 1* 

\A"A a + a^c" 



Table 1: Gauges and their renormalizable dimension two operator 



that explicit calculations outside the Landau gauge are almost prohibitive, as a localization proce- 
dure looks quite hopeless. 



Nevertheless, in several other gauges, we have shown that other dimension two, renormalizable, 
local operators exist. We generalized the LCO method and showed that these operators condense 
and give rise to a dynamical gluon mass, see Table 1 and [34, 35, |q]. In the maximal Abelian 



gauge, it was found that only the off-diagonal gluons A^ acquire a dynamical mass, a fact qualita- 



tively consistent with the lattice results from [|37 



We have been able to make some connection between the various gauges and their dimension two 



operators by constructing renormalizable interpolating gauges and operators [34, 39]. However, 



these operators are explicitly gauge dependent, hence there does not seem to exist a clear relation 
with the gauge invariant operator A^ in . 



Recently, the issue of Gribov copies has also been addressed in the maximal Abelian gauge [40, H]]. 



5. Another gauge invariant dimension two operator 

Recently, we have considered the fact that, perhaps, another gauge invariant dimension two 
operator might be of some relevance. We took a look at 



ab 



(5.1) 



This operator was already introduced by Jackiw and Pi during their analysis of a dynamical mass 



generation in 3-dimensional gauge theories [ }42[ ] . 



We can add the operator (|5JJ) to the Yang-Mills action as a mass term via 



with 



m 



id 



(5.2) 



(5.3) 



As we have discussed in [}43[], the action ( |5.2[ ) can be localized by introducing a pair of complex 
bosonic antisymmetric tensor fields, (fi^ v ,fi^ v ), and a pair of complex anticommuting antisym- 
metric tensor fields, (G^ V )^v)> belonging to the adjoint representation, according to which 



-Se 



DBDBDGDG exp 



d A xB a iiV Lf^D h ^B c iiV 
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Doing so, we obtain a classical local action which reads 

Sym + Sbg + S m , 

where 

SbG = 4 / ^ % { B ^v Da a D a BC nv ~ G ^v D ct &D ct G li 
Sm = l — fd 4 x(B-B) a .F" 



(5.4) 



(5.5) 



4 J " - \~ "/(iv'/iv 
which is left invariant by the gauge transformations 



(5.6) 



-Dfco b , 



SBl v = gf abc co% v , 8B; v =gf bc co% v , 
SG£ V = gf abc (o b G^ v , 8G; v =gf abc (0% v . 



(5.7) 



In order to discuss the renormalizability of (p.5|), we relied on the method introduced by Zwanziger 
in [22, 23] to discuss the renormalizability of the nonlocal horizon function ( 3.10p . Instead of 



using ( p.5| ) with m coupled to the composite operators B C ^ V F^ V and B^ V F^ V , we introduce 2 suitable 
external sources V pailv and V pa ^ v and replace S,„ by 



1 



^ d x [VapuvBopF^y V GpnvB% p F^ v j . 

At the end, the sources V aptlv {x), V aptlv {x) are required to attain their physical value, namely 

im 

phys ^ phys 2 



(5.8) 



V 



OpjlV 



{8 a ^8p V 8 av 8piij , 



(5.9) 



so that (5.8) reduces to S m in the physical limit. 



We assume the linear covariant gauge fixing, implemented through 

Sgf = J d 4 x (^b a b a + b a d^ +^d^D ab c 
In p3[], we wrote down a list of symmetries enjoyed by the action 

Sym + Sbg + Sgf , 



(5.10) 



(5.11) 



i.e. in absence of the sources. Let us only mention here the BRST symmetry, generated by the 
nilpotent transformation s given by 



sAl = -Dfc b , sc a = a -f bc c a c b , 

sB% = gf^cXv + G^ , sB; v =gf b ^c% v , 
= gf abc c b G^ v , sG; v =gf abc c% v +B% , 
sc" = b a , sb a = 0, s 2 =0. 



(5.12) 
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It turns out that one can introduce all the necessary external sources in a way consistent with 
the starting symmetries. This allows to write down several Ward identities by which the most 



general counterterm is restricted using the algebraic renormalization formalism [|43|]. After a very 

(5.13) 



cumbersome analysis, it turns out that the action ( p.5| ) must be modified to 



with 



S c i 



d 4 x 



1 f^a 



- -m 2 X\ (B^ v B^ v -G^ v G^v) + ' n2 ^f ( B % ~ b 1iv) 



+ 



A 



abed 



16 



H R 

D fiv D p.v 



32 



(5.14) 



in order to have renormalizability to all orders of perturbation theory. We notice that we had to 
introduce a new invariant quartic tensor coupling X abcd , subject to the generalized Jacobi identity 



jvnan j^mbed _|_ jmbn yarned _|_ jrmcn^abmd _|_ jrmdn^ 



rmdn i abem 



0. 



(5.15) 



and the symmetry constraints 



j^abed j^cdab 

j^abed ^bacd 



(5.16) 



as well as two new mass couplings X\ and A3. Without the new couplings, i.e. when X\ = 0, 
A3 = 0, X ahcd = 0, the previous action would not be renormalizable. We refer to [43, 44] for all 



the details. We also notice that the novel fields B a ^ v , B^ v , G^ v and G^ v are no longer appearing at 
most quadratically. As it should be expected, the classical action S c i is still gauge invariant w.r.t. 



The BRST transformation (5.12) no longer generates a symmetry of the gauge fixed action S p h ys . 
However, we are able to define a natural generalization of the usual BRST symmetry that does 



constitute an invariance of the gauge fixed action ( |5.13| ). Indeed, after inspection, one shall find 
that 



with 



sS 



phys 








(5.17) 



-Dfc b 



sc 



lfabc c a c b 



gf abc c% v , 7B a uv =gf abc c% v , 



gf bc c b G 



lb a =0. 



7G 



gf abc c% v 



(5.18) 
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In 44], we also calculated explicitly various renormalization group equations to two loop 
order, confirming the renormalizability at the practical level. Various consistency checks are at our 
disposal in order to establish the reliability of these results, e.g. the gauge parameter independence 
of the anomalous dimension of gauge invariant quantities or the equality of others, in accordance 



with the output of the Ward identities in [}43[]. Furthermore, we proved in [ ]44[ ] the equivalence of 
the model ( 5.13 ) with the ordinary Yang-Mills theory in the case that m = 0, making use of the 
nilpotent transformation 









= B % 


8 S (rest) 


= o, 



5 S G« V = 0, 



(5.19) 



which generates a "supersymmetry" of the action S™^. 
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